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A search for sidereal variations in the non-Newtonian force between two tungsten plates separated 
at millimeter ranges sets experimental limits on Lorentz invariance violation involving quadratic 
couplings of Riemann curvature. We show that the Lorentz invariance violation force between two 
finite fiat plates is dominated by the edge effects, which includes a suppression effect leading to lower 
limits than previous rough estimates. From this search, we determine the current best constraints 
of the Lorentz invariance violating coefficients at a level of ICE 8 m 2 . 
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Local Lorentz invariance is at the foundation of both 
the Standard Model of particles physics and General Rel¬ 
ativity (GR), however, the later theory is formulated as 
a classical theory, which demands some changes in its 
foundational structure to merge gravity with quantum 
mechanics. Even if local Lorentz invariance is exact in 
the underlying theory of quantum gravity, spontaneous 
breaking of this symmetry may occur, leading to tiny ob¬ 
servable effects [l], [2|. On the other hand, Lorentz viola¬ 
tions could also be large but "hard-to-see", despite many 
experiments to date setting very tight bounds across 
many physical sectors [H. This would occur if the Lorentz 
invariance violation is "countershaded" as pointed out in 
ref [U ■ Thus in general, the investigation of local Lorentz 
invariance violations in the spacetime theory of gravity 
is a valuable tool to probe the foundations of GR|5, b| 
without preconceived notions of the numeric sensitivity. 

Recently, the studies of Lorentz invariance violation 
in the pure-gravity sector shows that general quadratic 
curvature coupling will lead to interesting new effects in 
short-range experiments that could have escaped detec¬ 
tion in conventional studies to date Q- Accordingly a 
crude estimation has been made on the possible con¬ 
straints for these types of Lorentz invariance violation^, 
predicted to be tested to the level of ICR 8 to ICE 10 m 2 us¬ 
ing short-range experiments, such as the EotWash 9], 
Wuhan [Tol.lTlJ|. and Bloomington [T^| experiments. How¬ 
ever, for the latter two results, edge effects are not con¬ 
sidered properly. In this work we obtain the best cur¬ 
rent constraint of the Lorentz invariance violation at level 
10 -8 m 2 from a detailed reanalysis of prior data taken 
from the Wuhan experiment (HUST-2011). 

Effective field theory is a powerful and unique tool for 
investigating physics at attainable scales, and is suited 
for exploration of local Lorentz invariance in gravity. For 
centuries after Newton’s Principia, our experimental un¬ 
derstanding of gravity remains in some respects rernark- 
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ably limited. We are confident that Newton’s law de¬ 
scribes the dominant physics in long-range gravity and 
GR provides accurate relativistic corrections. However, 
in short-range gravity, it is presently unknown whether 
gravity obeys Newton’s law, and many models attempt¬ 
ing to unify gravity and the other fundamental forces in 
the same theoretical framework predict modifications of 
gravity, i.e., a deviation of the Newtonian 1/r 2 law. The 
experiments in short-range gravity are well suited for ex¬ 
ploration of local Lorentz invariance for spacetime-based 
gravitation. Here we use the data from a shot-range ex¬ 
periment in HUST-2011 to perform a search for short- 
range Lorentz invariance violation involving quadratic 
couplings of Riemann curvature. 

A quantitative description of Lorentz invariance viola¬ 
tion in the pure-gravity sector with quadratic couplings 
of Riemann curvature is given by the Standard-Model 
Extension (SME) Q. Using effective field techniques, 
the coupling leads to perturbative corrections to Newton 
gravity, which are inverse quartic and vary with orienta¬ 
tion and time. The corresponding perturbative correc¬ 
tions between two test masses mi and m 2 in the general 
quadratic curvature coupling of the SME is given by 


V lv (x)=-G- 
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where x = X\ — x 2 is the vector separating mi and m 2 , 
x 3 is the projection of the unit vector along x in the jth 
direction. {k e s)jkim is a set of 81 coefficients for Lorentz 
invariance violation in standard laboratory, which is to¬ 
tally symmetric with indices j, k , l , m ranging over the 
three spatial direction, representing 15 independent ob¬ 
servables for Lorentz invariance violation. These coeffi¬ 
cients can be taken as constant on the scale of the solar 
system. Therefore the canonical frame adopted for re¬ 
porting results from experimental searches for Lorentz in¬ 
variance violation is the Sun-centered frame jT3 UH , with 
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Z axis along direction of the Earth’s rotation and X axis 
pointing towards the vernal equinox. The zero in time 
is taken to be the vernal equinox at March 20, 2000 at 
7:35 a.m. Universal Time. Neglecting the Earth’s boost, 
which is older 10” 4 , the transformation from the Sun- 
centered frame ( X , Y. Z ) to the laboratory frame (x, y, z) 
involves a time-dependent rotation R? J . Taking the lab¬ 
oratory z axis pointing to the local zenith and the angle 
between x axis and local south being 9, the rotation ma¬ 
trix is 


R jJ 


cos 9 sin 9 0 
— sin 9 cos 9 0 


0 0 1 

cos x cos cj©T cos x sin w©T — sin x 
— sinw©T cosw©T 0 
sinxcosw©T sinxsinw©T sinx 


( 2 ) 


where x is the colatitude angle of the laboratory and 
w© ~ 27r/(23h56min) is the Earth’s sidereal frequency. 
The angle 9 is constant depending on the orientation 
of the apparatus in the experiment (For the experiment 
in Wuhan, we take 9 = —ir/2.). The sidereal time 
T-dependent coefficients ( k e g)jki m are thus related to 
the constant coefficients {k e &)j klm in the Sun-centered 
frame by 

(kes)jklm = R jJ R kK R lL R mM (hit)jKLM (3) 


Therefore, the inverse-cube potential for Lorentz invari¬ 
ance violation in Eq.(l) is oscillatory with T and includes 
components up to the fourth harmonic of w©, which leads 
to striking signals in short-range experiments. Equiv¬ 
alently, the Lorentz invariance violation force between 
two plates can be expected to vary with frequencies up 
to and including the fourth harmonic of w©. 

Most Inverse-Square Law tests use planar geometry to 
search for the non-Newton gravity. Noting that the New¬ 
tonian force at any point above an infinite plane of uni¬ 
form mass density is constant, the planar geometry can 
be effective in suppressing the Newtonian background rel¬ 
ative to the putative short-range effect. However, it also 
suppresses the Lorentz invariance violation signal pro¬ 
duced by Eq. (1). For instance, the Lorentz invariance 
violation force between a flat plate and an infinite plane is 
zero (shown in Fig. 1) due to the fact that the sum of the 
three coefficients in Eq.(l) is zero. The shape and edge 
effects play an important role in determining the sensi¬ 
tivity of the experiment to the coefficients for Lorentz 
invariance violation and is explained in following. 

Let’s consider the Lorentz invariance violation interac¬ 
tion in Eq.(l) between a flat plate of area Ap, thickness 
tp, and density pp at distance d from an infinite plane of 
thickness t and density p. Assuming the normal vector 
of the plates is in y axis direction as shown in Fig. 1, 
the Lorentz invariance violation acceleration on a point 
above the infinite plate with distance d is 


a v LV ( d ) = Gp 
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where the integration of the direction dependence 
k{9 p , ipp) can be carried out over the entire volume of 
the infinite plate. The integral in Eq(j4j) , shows that the 
Lorentz invariance violation force between the plates is 
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FIG. 1: The inverse-cube potential between two plates. 
Lorentz invariance violation force between a flat plate of area 
Ap, thickness tp, and density pp at distance d from an infi¬ 
nite plane of thickness t and density p is zero, which means 
the Lorentz invariance violation force between two finite flat 
plates is dominated by the edge effects. 


For the finite plates used in most experiments, one 
must use numerical integration. Noting that the mea¬ 
sured signal in experiments is usually the variation of the 
force when the gap between two plates changes from d m i n 
to d max , the variation of the Lorentz invariance violation 
signal can be roughly estimated as 


A Fly = Fl v (d min ) - F[ v (d max ) ~ eAC(k e s)jkim (6) 
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and the dimensionless parameter e represents edge ef¬ 
fects, typically of order 10 -2 or ~d/ y/~A in most inverse 
square law experiments. 

We then consider the experiment for HUST-2011. The 
design and operation of the experiment are described in 
Ref. Here we summarize briefly the basic features 

as shown in Fig.2. Two pure tungsten plates with the 
same dimensions of 16mm x 16mm x 1.8mm used as the 
test mass (Wl) and the counterbalance one (W2), were 
adhered to the two sides of the I-shaped pendulum sym¬ 
metrically. The arm length of the test mass is 42mm. The 
source mass platform, facing to the pendulum, was also 
designed as I-shaped structure. A larger tungsten source 
mass (W3) with dimension 20.8mmx 20.8mmx 1.8mm 
was adhered to the one side of the I-shaped structure, 
opposite to the test mass Wl. The non-Newtonian force 
to be measured is between Wl and W3 for separations 
ranging from 0.4 to 1.0 mm, where the net torque change 
of the Newtonian force was counteracted by a thicker 
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FIG. 2: Schematic drawing of the experiment HUST-2011. 
The pendulum twist was measured by an autocollimator, and 
controlled by applying differential voltages to the two capaci¬ 
tive actuators (not shown here). The separation between the 
counterbalance masses W2 and W4 (« 4mm) is larger than 
that between the test mass Wl and source mass W3. The 
separation between Wl and W3 was modulated by driving 
a motor translation stage with frequency f s , and the sensi¬ 
tivity of the pendulum was calibrated by rotating the copper 
cylinder at frequency f c . 


tungsten counterweight mass (W4) with dimensions of 
16mmx 16mmx7.6mm, adhered on the other side of the 
source mass platform, located behind W3 about 4 mm 
in y axis direction. The normal vectors of the plates are 
all in y axis direction in local laboratory frame (Fig.l 
and Fig.2). Taking Ap=16mmx 16mm, tp=t=1.8mm, 
arm length l 42mm, d m i n =0.4mm, and d max =1.0mm, 
we obtain AChust^ ~ 9-42 x 10 _7 Nm/m 2 according 
to Eq.(7). The numerical integration gives the torque 
for Lorentz invariance violation between test mass and 
source mass, and it is determined to be 


At 


LV 


-0.06(fc eff ) rara + 0.16 (k eS ) XX yy + ■■■ 

xlO 7 Nm/m 2 << ACnusT^eff) 
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which means the edge effect suppression factor is approx¬ 
imately e ~0.01. Taking the sensitivity of the torque as 
0.8xl0~ 16 Nm, the constraint of ( k e g ) can be roughly 
estimated at level of 0.5 xlO -8 m 2 . Then the detailed 
analysis of sidereal variations of torque can set the ex¬ 
perimental limits on Lorentz invariance violation param¬ 
eters, (fceff)jifLM, in the Sun-centered frame. 

In the experiment for HUST-2011, the motion of 
the pendulum was controlled by using a proportional- 
integral-differential feedback system. The separation be¬ 
tween the tungsten test and source masses was modu¬ 
lated by driving a motor translation stage with frequency 
f a = 27t/500s, and the sensitivity of the pendulum was 
calibrated by rotating a copper cylinder at frequency 
f c = 27t/400s. The recorded data involves 12 sets, ac¬ 
quired from August to October in 2011. Each set has a 
duration time of approximately one day. To extract the 
expected signal, each data run was broken into separate 
“cuts” containing exactly the lowest common multiple of 
the dual-modulation periods, AT =2000s. For each “cut”, 


and then the non-Newton torque was fitted with 


T measured( T ) = TLv(T) COs(2t T f s T + ip) (9) 

with the initial phase p for the separation modulation, 
which can be determined according to the operation pro¬ 
cedure in the experiment. For averaging times of AT 
=2000s, the sidereal variation of tlv(T) can be neglected 
since w ffi AT << 1. Thus, we take rpy (T) as constant ap¬ 
proximately in each “cut”. Fig.3 represents the sidereal 
variation of TLv(t) varying with time t. Each point rep¬ 
resents the mean of a 2000 second measurement without 
showing the error bar, which is dominated by statistical 
uncertainty. In the Sun-centered frame, the time scale is 
reset by t = 0 when the Sun passes the vernal equinox on 
2011, March 20, 23 h 21 min universal time (UT). In the 
laboratory frame, the time scale is set by T = 0 when the 
east direction coincides with the Y axis. The difference 
between the two time scales can be written as a phase 
difference w®(T — t) ~ —1.35 for this experiment (east 
longitude 114.4° and x = 59.5°). 

According to Eqs. (2) and (3), tlv(T) can be further 
expressed as a Fourier series in the sidereal time T 
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tlv(T ) = c 0 + ^2 i c m cos(mw©T) + s m sin(mw©T)] 

m=1 

sin(mcu ffi AT/2) 

mw® AT/2 ’ 1 1 

in which we have considered corrections due to data 
averaging. For the sinusoidal signal with angular fre¬ 
quency men®, the data average over AT intervals will 
lead to the attenuation of amplitude by a small amount 
1 — sln r ^^r/' 2 2 ^ ’ w hich i s approximately zero. The 
Fourier amplitudes in this expression are functions of the 
coefficients (fc e ff) jklm, the test mass and source mass 
geometry, and the laboratory colatitude. For example, 
the coefficients cq and Ci related to (k e s) jklm can be 
obtained by numerical integration, 


co = [— Q.Ql2{k e s)xxxx — 0.012(/c b s)yyyy 
+0.025 (k e s)zzzz — 0.024(/c ef r)A-jfyr 

+0.0ll(k e fi)xXZZ + 0.011 \k e s)YYZz] 

xl0 _7 Nm/m 2 , . 

Ci = [— 0.085(fc e ff)AA'AZ “ 0.002 (k e ff)YYYZ 
—0.006(_fc e ff)zxzA' — 0-Q32(k e fi)zzzY 
—0.002(k e g)xxYZ — 0.085(/c e ff)rYA'z] 
x 10” 7 Nm/m 2 . 

The systematic error of T£y(T) is about 0.9xl0~ 16 
Nm, which comes from the dimensions and positions of 
the test and source mass, as well as the torque calibra¬ 
tion. The largest part of systematic error is due to the 
uncertainty of Newtonian torque. However, it is indepen¬ 
dent on sidereal time T, and only affects the coefficient co- 
Therefore, the uncertainty of co is dominated by the sys¬ 
tematic error, and the uncertainty of Ci and s, (i=l,2,3,4) 
is dominated by the statistical error. The nine compo¬ 
nents ci and Si extracted from the experiments by Eq.(10) 
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FIG. 3: Data for TLv(t) from the HUST-2011 short-range 
experiment. 


TABLE I: Fourier amplitudes in 10 16 Nm units. 


Mode 

Numerical value 

co 

-0.22±0.95 

Cl 

0.13±0.22 

Sl 

-0.40i0.23 

C2 

-0.04i0.22 

S2 

0.20i0.22 

C 3 

-0.30i0.22 

S3 

0.25i0.23 

C4 

-0.06i0.23 

S4 

0.05i0.23 


are listed in the Table I. Each of them is functions of the 
coefficients ( k e s)jKLM > such as in Eq.(ll). 

However, the nine signal components are insufficient to 
constrain independently each of the 15 degrees of freedom 
in (kes)jKLM- Following standard practice in Ref. [12), 
we can obtain maximum-sensitivity constraints on each 
component of (fc e ff)jifiM in turn by setting the other 14 
degrees of freedom to be zero. The resulting measure¬ 
ments and 2er errors on each independent component of 
(kes)jKLM are displayed in Table II, which is in agree¬ 
ment with the rough estimation at level of 0.5 xlO -8 m 2 
in Eq.(8). 

The Table II represents the best measurements of pos¬ 
sible short-range Lorentz invariance violation involving 
perturbative inverse-quartic corrections to Newton’s law 
and hence of possible quadratic curvature couplings vio¬ 
lating local Lorentz invariance. As seen above, an impor¬ 
tant feature of short-range tests of local Lorentz invari¬ 
ance violation in gravity is that the sensitivity to Lorentz 
invariance violation is dominated by edge effects. As a 
comparison let us consider the experiment of Indiana- 
2012 [12]. In the experiment, each of the two masses 
is a planar tungsten oscillator of approximate thickness 
250 /rin, separated by a gap of about 80 /./m Ena. 
The force-sensitive ‘detector’ mass is driven by the force- 


TABLE II: Coefficient values (2a) from the experiment 
HUST-2011 


Coefficient 

Value (10~ 8 m 2 ) 

(k e s)xxxx 

-0.2i2.8 

( k e ff)yYYY 

0.4i2.8 

(kefl)zzzz 

-0.9i7.7 

(keff)xXXY 

0.4il.3 

(keff)xxxz 

-0.1i0.5 

\keSJYYY X 

0.6il.3 

(kefl)YYYZ 

0.4i0.5 

(k e g)zzzx 

-1.3il.4 

( k e ff)zZZY 

-0.2il.3 

(k e ff)xXYY 

-0.1il.7 

(k e ff)xXZZ 

-0.2il.0 

{k e s)YYZZ 

0.2il.0 

(k e fl)xXYZ 

0.5i0.5 

(k e s)YYXZ 

-0.2i0.5 

{k e s)zZXY 

-0.2i0.5 


generating ‘source’ mass at a resonance near 1 kHz. De¬ 
tector oscillations are read out via a capacitive trans¬ 
ducer probe coupled to a sensitive differential amplifier, 
with the signal fed to a lock-in amplifier referenced by 
the same waveform used to drive the source mass. The 
overlapping area of the test mass and source mass is 
j4 p =5.7mmx5.08mm, separated by a gap 60 fim. for G? m ; n 
and 100 (im for d max . According to Eq.(7), we obtain 
ACindiana ~ 15-1 x 10~ 7 N/m 2 . The preliminary numeri¬ 
cal integration gives 


A F, 


LV 


O.l(fceff) 

xxxx + 0 - 5 (k e ff) XX yy + ■ ■ • 

X10 N/m << AClnUiana^effi 


( 12 ) 


which means the edge effect suppression factor is approx- 
imatively e ~ 0.03. Therefore, the sensitivity of the ap¬ 
paratus to the coefficients (fc e ff) jklm can be roughly es¬ 
timated as the ratio of the thermal-noise force (~10 fN) 
to the scale (~0.05 pN/m 2 ) in Eq.(12), resulting in a 
level of 2xl0~ 7 m 2 . 

In conclusion, in short-range tests of local Lorentz in¬ 
variance violation in the gravity sector, Lorentz invari¬ 
ance violation force between two finite fiat plates is dom¬ 
inated by edge effects, which adds an extra suppression 
factor. Nevertheless, this work presents the best cur¬ 
rent constraints of Lorentz invariance violation involving 
quadratic couplings of Riemann curvature in the pure 
gravity sector at level 10” 8 m 2 . 
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